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CO ' Abstract 

<N : 

We study photon-induced coherent production of T in ultraperipheral heavy ion 

collisions at LHC and demonstrate that the counting rates will be sufficient to measure 

qq | nuclear shadowing of generalized gluon distributions. This will establish the transition 

7— I ■ from the regime of color transparency to the regime of perturbative color opacity in 

an unambiguous way. We argue that such measurements will provide the possibility 

to investigate the interaction of ultra-small color dipoles with nuclei in QCD at large 

C*~> ■ energies, which are beyond the reach of the electron-nucleon (nucleus) colliders, and 

will unambiguously discriminate between the leading twist and higher twist scenarios 

*^ of gluon nuclear shadowing. 

P-i! 

Jh ; 1 Introduction 

> ■ 

High energy coherent vector meson photo- and electroproduction provides new tools for 
understanding of many exciting phenomena in the QCD physics of the hadron-nucleus in- 
teractions. The photon wave function comprises a variety of hadronic components as well 
as a direct contribution of the quark-antiquark pairs. In the high energy regime, the photon 
transforms into a qq pair long before the nuclear target. The transverse size of the produced 
qq pair is controlled by the mass of the quark (r oc m" 1 , provided m q ^> Aqcd)- Corre- 
spondingly, the scale of virtuality in the considered processes is Q 2 > vr? . The strength of 
the high-energy interaction of small dipoles in QCD depends on the size of the area occu- 
pied by the color field within the interacting objects: the smaller the size, the weaker the 
interaction P]I2]- Probably the most sensitive indicator of the size of the interacting objects 
is the t-dependence of exclusive vector meson electro- and photoproduction. The current 
HERA data on the vector meson electroproduction on the proton target, see e.g. |3 a , are 
consistent with the prediction of jlj that the t-slopes of the p and J/ip production production 
cross sections should converge to the same value with the rate of convergence consistent with 
the estimate of This indicates that at small x, x — [My + Q 2 )/s (x is related to the 
momentum fractions of the exchanged gluons - see the discussion of Eq. (Jlj; s is the invari- 
ant energy for 7 — N scattering), the configurations of much smaller size than the average 
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light meson size dominate J/ip production for all Q 2 and p production for Q 2 > 5 GeV 2 . 
In the kinematics of small x and large Q 2 , the QCD factorization theorem predicts that 
the cross sections of coherent production of vector mesons off nuclei are proportional to the 
square of the nuclear gluon parton density 1 Ga(x, Q 2 ) [UEUE]. It is reasonably well known 
from various data analyses that Ga(x, Q 2 ) is not shadowed and even may be enhanced at 
0.02-^0.03 < x < 0.2: G A {x,Q 2 ) > AG N (x,Q 2 ). One should note that the smalkr behavior 
of nuclear parton distributions is crucial for various aspects of small-x dynamics with nuclei. 
The nuclear parton distributions have been analyzed both phenomenologically, by fitting 
to the available data (HI E| (which for Q 2 > 2 GeV 2 are available only for x > 10~ 2 ), and 
theoretically, by using various models, see Ref. (THj for the review. These analyses suggest a 
rather significant nuclear shadowing assuming that the leading twist effects dominate. How- 
ever, the F 2 a/F 2 d data are restricted to rather moderate Q 2 , where higher twist effects due 
to e.g. diffractive vector meson production are not negligible. 

On the theoretical end, there are two types of models. One class focuses on sufficiently 
large virtualities, where leading twist dominates. It assumes that the soft physics (aligned 
jet model) which dominates at the starting scale of evolution for quark shadowing should 
be present in the gluon sector as well, so that the gluon shadowing effects at the initial 
evolution scale of a few GeV 2 should be comparable to that of quarks JTJ [H] . More recently 
a more quantitative approach became possible. It is based on the combination of the Gribov 
theorem connecting diffraction and shadowing [I2j, arid the Collins factorization theorem 
for hard diffraction in DIS [13]. As a result, this approach effectively takes into account 
the nonconservation of the number of particles in QCD evolution and the restoration of the 
gluon fields by the small color singlets. When combined with the HERA hard diffractive 
data and quasieikonal modeling of the multiple rescattering terms, it leads to the prediction 
that the gluon shadowing is significantly larger than that for quarks at the starting scale of 
QCD evolution of Q 2 ~ 4 GeV 2 . 

Another group of the models J3] aims to model the effects of the modification of nuclear 
gluon field both in the regime of high gluon densities, where the decomposition over the 
twists is not applicable and at large Q 2 where the leading twist dominates. It assumes that 
parton densities are saturated at Q 2 < Q 2 (Q 2 (x) is the saturation scale). An important 
feature of these models is that they do not lead to the leading twist shadowing at Q 2 ^> Q 2 , 
see e.g. discussion in [Tallin]- in particular, in these models the effect of nuclear shadowing 
for the interaction of a small dipole is modeled in the frozen impact parameter eikonal 
approximation with the dipole coupled to the unscreened nucleon gluon field, leading to the 
shadowing oc xGn(x, Q 2 )/Q 2 - 

A large leading twist gluon shadowing would significantly slow down the increase of the 
cross section of the small dipole-nucleus interaction with energy (which is proportional to 
the gluon density) and, hence, would significantly extend the x range, where unitarity for 
the interaction of small dipoles with nuclei is not violated in the DGLAP approximation. As 
a result, this may strongly affect the pattern of the possible onset of the regime of the black 
body limit (BBL). 

In the arrange, where nuclear gluons are not shadowed, one expects the regime of color 

1 Note that the corrections due to the skewedness effects at large Q 2 and small x are calculable in terms 
of the QCD evolution equation for the generalized parton distributions. 
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transparency for coherent production of onium states. A leading twist shadowing would 
lead to the onset of the perturbative color opacity with the decrease of x. However, this 
shadowing merely slows down the increase of the gluon density with the decrease of x. 
Hence, at sufficiently small x, the perturbative color opacity regime will be also violated. 
This may lead to the onset of BBL or another nonlinear regime of QCD. 

The crucial questions for understanding the interplay among the above mentioned phe- 
nomena are the following. What is the kinematical region where the squeezing of the qq 
dipole takes place, which justifies the application of perturbative QCD? What are the lead- 
ing mechanisms of nuclear shadowing at very small x, where the small size dipole-nucleon 
interaction become large? Can the gluon saturation effect at small x stop the increase of the 
cross section before reaching the BBL. 

In order to establish the effect of nuclear opacity and clarify the role of the leading and 
higher twist mechanisms of nuclear shadowing in an unambiguous way, it is necessary to 
investigate the interaction of a small dipole with the nuclear medium and check whether the 
amplitude of the interaction is reduced as compared to the color transparency expectation of 
the amplitude being proportional to the atomic number. Production of light vector mesons 
at small x in DIS by the longitudinally polarized photons and photoproduction of onium 
states seem to be optimal for these purposes. The studies of the dipole model in Ref. [T7] 
suggest that one is safely into the perturbative domain at x < only if the transverse 
distances are below d ~ 0.3 fm. For the light mesons this requires Q 2 > 15 GeV 2 and, 
hence, an eA collider with energies of HERA and large luminosity. The photoproduction 
of J /if) appears to be the borderline case - the amplitude seems to be getting a significant 
contribution from the transitional region between perturbative and nonperturbative QCD. 
Although J / ip photoproduction suits well the studies of the onset of color opacity, in general 
it may be rather problematic to distinguish the perturbative regime, where the amplitude 
is proportional to the gluon density, from the regime, where the fast increase of the gluon 
density results in the violation of the leading twist approximation and the transition to the 
black body limit (BBL) takes place. In Ref. ^H] it was suggested to study the onset of BBL 
in the coherent photoproduction of dijets from heavy nuclei since in this case the role of 
the leading twist shadowing is negligible. On the contrary, the coherent T photoproduction 
off heavy nuclei seems to be essentially unique for the investigation of the onset of the 
perturbative color opacity in the leading twist (LT) shadowing regime. In this case, the 
interaction is dominated by such small transverse interquark distances d t , d t ~ 0.1 fm, that 
physics of interactions is definitely perturbative up to very large energies. For example, the 
cross section of the interaction of a dipole with d t ~ 0.1 fm is a(d t = 0.1 fm, x = 10~ 4 ) = 3.5 
mb and is even smaller for larger x. 

In order to probe the onset of perturbative color opacity in the T production, one would 
need to reach x <C 10~ 2 . This would require an electron- nucleus collider in the HERA 
energy range with a luminosity per nucleon comparable to the luminosity in the electron- 
proton collisions, which appears to be rather problematic for HERA. At the same time there 
are plans for studying the photoproduction off nuclei at LHC using ultraperipheral heavy 
ion collisions, including the studies of the coherent onium photoproduction. For the review 
and extended list of references, see Ref. [T5j . 

Recently it was demonstrated {20 that the yield of J / if) in the coherent production in the 
kinematics of the ultraperipheral ion collisions at LHC is suppressed at the central rapidities 
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by a factor of 6 compared to that in the impulse approximation. Such a significant effect of 
the suppression would clearly signal the revealing of the color opacity phenomena. However, 
as mentioned above, the gluon virtualities in this case are on the borderline between the 
perturbative and nonperturbative QCD regimes (Qlff ~ 3 4- 4 GeV 2 ) and the leading twist 
approximation might not be accurate enough for the evaluation of the scattering amplitude, 
especially since the nonperturbative region could probably give a significant correction |21j . 

In summary, T photoproduction stands out as the simplest probe of the propagation of 
small color dipoles, d t = 0.1 fm, at LHC energies, which will allow to investigate the nuclear 
gluon fields at the effective scale Q 2 ~ 40 GeV 2 . In this case, the estimates of nuclear 
shadowing via the leading twist mechanism based on perturbative QCD used in Ref. [20J, 
will be much better justified. The main objective of this work is to present the results of a 
study of the coherent production of T in the ultraperipheral collisions of heavy ions in the 
kinematics of LHC. 



2 Production of T in ultraperipheral collisions 

The cross section of the coherent T production integrated over the transverse momenta of 
the nucleus, which emitted a photon, can be written using the standard Weizsacker- Williams 
approximation [22] in the form: 

da(AA -> TAA) Ar . , , . Ar . , . . 

j- = N 7 (y)cr jA ->rA{y) + N 1 {-y)a 1 A-,rA{-y) ■ (1) 

Here y = i n 2 "t m NM v ^ g ^ e c m rapidity of the produced T; 7 is Lorentz factor; cr(y) is the 
cross section for the photoproduction of T; N 7 (y) is the flux of the equivalent photons 

Z 2 a f , Q ,„ /T *, 1 „„„,,„, 1 

1 



N(y) = ^J d 2 bT AA (b)±X 2 [K 2 (X) + -Kl(X)] , (2) 



where K (X) and K\{X) are modified Bessel functions of the argument X = bmyt y / (27); b 
is the transverse distance between the centers of the colliding nuclei. The standard Glauber 
profile factor, 

00 

T AA (b) = ex.p(-a NN J dz J d\p A (z,b 1 )p A (z,b-b 1 ) S J , (3) 

—00 

accounts for the inelastic strong interactions of the nuclei at impact parameters b < 2R A 
and, hence, suppresses the corresponding contribution to the T photoproduction. In our 
calculations we use the nuclear matter density p A (z, b) obtained from the mean field Hartree- 
Fock-Skyrme (HFS) model, which describes many global properties of nuclei [25 as weu as 
many single-particle nuclear structure characteristics extracted from the high energy A(e, etp) 
reactions [25] . 

The amplitude of the T photoproduction (necessary for the calculation of <j iA ^ta i n 
Eq. (^)) in the leading twist approximation is described by the series of the Feynman diagrams 
in Fig. [21 The QCD factorization theorem for exclusive meson photoproduction [3 El E] 
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Figure 1: High energy quarkonium photoproduction in the leading twist approximation. 



allows one to express the imaginary part of the forward amplitude for the production of a 
heavy vector meson by a photon, 7 + T — > V + T, through the convolution of the wave 
function of the meson at the zero transverse separation between the quark and antiquark, 
the hard interaction block and the generalized parton distribution (GPD) of the target, 
GT(xi,X2,Q 2 ,t min ), at t m i n ~ —x 2 m 2 N . The gluon light cone fractions Xi for the gluons 
attached to the quark loop satisfy the relation: 

x\ - x 2 = = x , (4) 

s 

where s = AE^u = AjurniN is the invariant energy for 7 — N scattering (En = ^rriN is 
the energy per nucleon in the cm. of the nucleus-nucleus collisions). If the quark Fermi 
motion and binding effects were negligible, x 2 C ii. Numerical estimates using realistic 
potential model wave functions indicate that for J/ip, x\ ~ \.5x,x 2 ~ x/2 [21], and that 
for T, x 2 /xi ~ 0.1 [26J. Modeling of the GPDs at moderate Q 2 suggests that, to a good 
approximation, Gt(x\,x 2 ) can be approximated by the gluon density at x — (x\ + x 2 )/2 
[3J 121]. F° r large Q 2 and small x, GPDs are dominated by the evolution from Xi(init) ^> X{. 
Since the evolution conserves X\ — x 2 , the effect of skewedness is determined primarily by 
the evolution from nearly diagonal distributions. 

Hence, we can approximate the amplitude of the T photoproduction off nucleus at k 2 = 0, 
presented by a series of diagrams in Fig. [21 as 

M( 7 + A - T + A) = M( 7 + N - T + A) ^ ( *' F A (t min ) , (5) 

ALr N {X, ty e ff) 

where Fa is nuclear form factor normalized so that Fa(0) = A; Q 2 j f (T) ~ 40 GeV 2 according 
to the estimates of jHj- It is worth noting that really we are not sensitive here to the precise 
value of Qlff since for large Q 2 , the gluon shadowing decreases slowly with Q 2 . Note that the 
skewedness effects, which are expected to be small at Q 2 ~ a few GeV 2 , become important 
in the T case. They appear to increase the cross section of the reaction 7 + p — > T + pby 
a factor ~ 2, see Refs. [2E1I2E]- Potentially, this could obscure the connection between the 
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perturbative color opacity effect and shadowing of the nuclear gluon densities. However, the 
analysis of |29j shows that the ratio of GPD on nucleus and on a nucleon for t = t min is a 
weak function of x 2 , slowly dropping from its diagonal value (xg = x\) with the decrease 
of X2- Overall this observation is in a agreement with the general trend mentioned above 
that it is more appropriate to do a comparison of the diagonal and non-diagonal cases at 
x — (xi + x 2 )/2. 

The effect of leading twist nuclear shadowing can be quantified by considering the dia- 
gram (b) in Fig. [2] All possible unitary cuts of the diagram lead to inelastic shadowing in 
hadron-nucleus total cross sections [12]. Applying these ideas to DIS on nuclei and using 
the factorization theorem for hard diffraction [HJ, coupled with the QCD analysis of the 
HERA data on diffractive DIS [3U], nuclear parton distributions can be predicted [T3] ITU]. 
In general, the theory of leading twist nuclear shadowing should be applied in two steps. 
Firstly, the nuclear gluon density distribution is calculated at the starting evolution scale of 
the effective momentum transfer Ql = 4 GeV 2 



(1 - in) 2 

G A (x, Ql) = AG N (x, Ql) - 87T^e 



1 + 7/ 2 



oo oo -^jp 

J d 2 b J dz\ J dz% J dxjp 



xg%{— , x P , Ql t mm )p A {b, Zl ) PA {b, Z2 ) e ^ m ^-^ e 



OO Z\ X 

\ a ef f(x,Ql)(l-iri) J dz PA (b,z) 



• (6) 



Here xjp, Ql, tmin) is the diffractive gluon density distribution of the nucleon; xjp is the 

Pomeron momentum fraction; n is the ratio of the real to imaginary parts of the elementary 
diffractive amplitude; x P = 0.03 is the cut-off parameter of the theory. As an input, we 
used the HI parameterization [30] of X]p, Qo, tmin)- The effect of the interactions with 

three and more nucleons of the target (graph (c) in Fig. [2]) is included in the attenuation 

factor T = exp(— a e jf(x, Qq)(1 — iT))/2 J dzp A (b, z)). The effective cross section a e ff(x,Ql) 

Zl 

accounts for the elastic rescattering of the produced diffractive state off the nuclear nucleon 
in the series of the rescattering diagram (graphs (b) and (c) in Fig. [2]) and is defined by the 
equation 



16tt 



a eff (x, Ql) - f dxjpg§( — , x P , Q 2 , t min ) . (7) 

[i + r] )Lr N {X,(^Q) J Xp 

Secondly, since the photoproduction of T corresponds to larger Q 2 scales ~ 40 GeV 2 , the 
gluon density distributions were evolved up to this scale using the NLO QCD evolution 
equations. In the limit of low nuclear density, when the effect of the attenuation factor T 
can be neglected, Eq. (JBJ) can be applied to evaluate G A at any scale Q 2 and not only at Qq. 
In the general case, Eq. (0) should be used only at rather moderate Ql and should serve 
as an input for QCD evolution: The modelling of multiple rescattering by the attenuation 
factor T ignores the fluctuations in the strength of the interaction around the average value 
u e ff, which is a good approximation only at moderate Qq. Indeed, at Ql ~ 3 4- 4 GeV 2 , the 
interactions are predominantly soft with small cross section fluctuations. The effect of cross 
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section fluctuations is automatically included in QCD evolution, which leads to violation of 
the Glauber - like structure of the expression for the shadowing for Q 2 > Q%. 

One can also address the question of the accuracy of the substitution of the ratio of the 
generalized gluon densities by the ratio of the diagonal parton densities at the normalization 
scale. In the case of scattering off two nucleons (only graph (b) in Fig. EJ), one can express 
the leading twist screening via the nondiagonal matrix element of the diffractive distribution 
function, g D (an analog of generalized PDF). It depends on the light-cone fraction, which 



the nucleon lost in the \in) and (out\ states: Xjp Xjp — X respectively, 
(xi — x)/ (xjp — x), and t, Q 2 . If we make a natural assumption that 



%l/ Pout 



g D {x 1 , x, xjp, Ql, t) = yj g D (f3 in} Q% x P , t)g D (f3 out} Q§, x P - x,t) , (8) 

we find that in the kinematics we discuss, the resulting skewedness effects are small nu- 
merically. They are smaller than the uncertainties in the input gluon diagonal diffractive 
PDFs. 

The discussion of the nuclear gluon distribution culminates in the fact that the cross 
section of the process — > TA can be readily obtained by squaring the amplitude in 
Eq. © 



T ~fA->TA(S) 



where —t 



(jcr 7 jv_, T jv(g, t min ) 
dt 



G, 



X ,Q 2 



effJ 



AG N (^,Q 2 eff ) 



2 U, 



dt 



d 2 bdze^\- iq ' z p A (b,z) 



(9) 



m 2 + 



(ft | 2 is the square of the vector meson transverse momentum and qi 



tunM^/s is the minimal longitudinal momentum transfer in the photoproduction vertex. 

The elementary cross section of forward T photoproduction in the energy range of interest 
is not known well experimentally. In our calculation, we used a simple parametrized form 



da. 



S,t) 



dt 



10~ A B r 



s 
so 



0.85 



(10) 



with the reference scale so = 6400 GeV 2 ; the slope parameter Br = 3.5 GeV~ 2 and the 
energy dependence, which follows from the calculations j2E| of the photoproduction of T in 
the leading log Q 2 approximation with an account for the skewedness of the partonic density 
distributions. The cross section is normalized so that the total cross section is in yh. 

The results of the calculations using Eq. (jH)) are presented in Fig. 13 The figure depicts 
the ratio of g 1 a-^ta calculated with Eq. (JHJ) to that calculated ignoring the effect of nuclear 
shadowing (setting Ga = AGn)- We find a rather strong suppression of the coherent T 
photoproduction cross section off nuclei. Note that UPC collisions at LHC are sensitive to 
the x-range of 3 x 10~ 4 < x < 5 x 10~ 3 . 

Since the significant effect of shadowing is far beyond the energies available in the fixed 
target experiments with the photon beams, the only opportunity to study this phenomenon is 
coherent production in UPC at LHC 2 . In this kinematics, the coherent events can be selected 

2 In principle, HERA can reach into the necessary kinematics. However, this would require both the 
approval of HERAIII program and development of the system necessary for high luminosity runs with heavy 
nuclei. 
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Figure 2: The shadowing effect in the T photoproduction off Pb and Ca. 
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Figure 3: The momentum transfer dependence in the coherent T production in the UPC at 
LHC The solid curve corresponds to the coherent cross section with accounting for nuclear 
shadowing; the dashed curve corresponds to incoherent T production. 
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Figure 4: The rapidity distribution for the coherent T production in Ca-Ca and Pb-Pb in 
UPC at LHC. The solid curve corresponds to the calculation including leading twist nuclear 
shadowing; the dotted curve corresponds to the calculation with the model of shadowing 
of Eskola et ai; the dot-dashed curve is the calculation in the eikonal dipole rescattering 
model; the dashed curve corresponds to the impulse approximation. 



10 



Approximation 


CaCa at LHC(7 = 3500) 


PbPb at LHC(7 = 2760) 


Impulse 


1.8 fib 


133 fib 


Glauber+Leading twist shadowing 


1.2 fib 


78 fib 



Table 1: Total cross sections of T production in UPC at LHC. 

by the requirement of anticoincidence with the signal in the Zero Degree Calorimeter( the 
requirement is that there are no neutrons in the final state) and by selecting the produced 
T with the small transverse momentum, which strongly suppresses the contribution of the 
incoherent diffraction. The specific transverse momentum distributions for both coherent 
and incoherent production are shown in Fig. El The cross section of the incoherent T 
photoproduction was estimated in the impulse approximation as 

da~fA^rx(s, t) _ ^ da^ N ^y N (s, t) 
dt dt 

This constitutes the upper limit for incoherent contribution. 

The rapidity distributions for coherent T production in the UPC with Ca and Pb beams 
are shown in Fig. 0] and the corresponding total cross sections are given in Table H For 
comparison in Fig. HI we also present results (the dotted curve) of the calculation using the 
model of Eskola et al. [Hj for the gluon shadowing, which predicts somewhat smaller suppre- 
sion of the T yield. In order to illustrate that that eikonal dipole models of rescatterings 
give a much smaller suppression, we present also the result of the eikonal calculation (the 
dot-dashed curve) using the analysis of the value of a(d, x) from jT7j(the result of [HT| for 
d ~ 0.1 fm are very similar, see comparison in |17j). It appears that a detector with good 
acceptance for T production for the currently discussed luminosities for heavy ion runs of 
LHC will collect enough statistics to measure the cross section of the discussed process with 
good precision. As seen from comparison to the calculations in the Impulse Approximation 
(Fig. HJ), the yield of T is expected to be suppressed by a factor of two at central rapidities 
due to leading twist shadowing. Practically, the only opportunity to draw the conclusion 
about the onset of the color opacity from the measurement of the T yield in the ultraperiph- 
eral ion collisions at LHC, which we see at the moment, is the comparison of the calculations 
with the data. Obviously until the rather large uncertainty in the absolute value of the 
7 + N — > T + iV cross section exists, such a conclusion cannot be considered as decisive. 
However, in the coming couple of years it is expected that this elementary cross section will 
be measured at HERA with much better accuracy. Besides, the energy dependence of the 
elementary cross section can be considered as well established theoretically. Hence, a mea- 
surement of the T yield at the edge of the rapidity distribution, where the effect of shadowing 
is still insignificant, could help to fix the uncertainty and verify the calculation of the cross 
section in the impulse approximation in all ranges of rapidities. Therefore, a comparison 
with the data at the central rapidities would provides us with an estimate of the shadowing 
effect. Note that there exists a procedure [223 t° separate the production of mesons at small 
and large impact parameters through the study of the nucleus break up, which would allow 
to determine the ratio of the cross section of T production by left and right moving photons 
for a given rapidity. Another possibility will be to measure the elementary cross section 
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at LHC in UPC collisions of protons with nuclei. Though this cross section is a factor of 
oc A 4 / 3 smaller, the expected luminosity in pA collisions is a factor of 100 larger than in AA 
collisions. 



3 Conclusion 

We found that the UPC program at LHC will provide a practical way to search for the onset 
of the perturbative color opacity via the study of the coherent T production at LHC and 
will decisively discriminate between the leading twist and higher twist scenarios of gluon 
shadowing at large Q 2 . 

This work was supported by GIF, Sofia Kovalevskaya Program of the Alexander von 
Humboldt Foundation and DOE. 
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